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Computer System and Method for Calculatine a Pharmacokinetic Behavior of a 
Chemical Substance in Insects 

The invention relates to a computer system for calculating a 
pharmacokinetic behavior of a chemical substarice in insects on the basis of a 
physiological model, as well as to a corresponding method and computer program 
5 product. 

Background of the invention 

Many physiologically based models for manmials are known from the 
prior art (for example, Chamick et al., J. Pharmacokin Biopharm. 23, 217 (1995)). A 
physiologically based model for caterpillars is furthermore known from the prior art 
10 (Greenwood et al, Pestic. Sci. 30, 97 (1990)). The insect is in this case described by 
compartments, which each represent an individual organ of the insect. The 
"interconnection" of the individual compartments is obtained from the known 
physiology of caterpillars. The essential parameters of this model are the rate 
coefficients for inter-compartmental mass transport, which determine the speed of 
15 the distribution, and the organ distribution coefficients, which give the concentration 
ratio in thermodynamic equilibrium between the respective organ and the 
haemolymph, which corresponds to the blood fluid in mammals. 1 

it is known from the prior art to use such compartmental models in order 
to describe an experimentally determined pharmacokinetic profile of a substance 
20 retrospectively, by matching the rate coefficients and the distribution coefficients (for 
example, Lagadic et al., Pestic. Biochem. Physiol. 45, 105 (1993) & Pestic. Biochem. 
Physiol. 48, 173 (1994)). 

It is therefore an object of the invention to provide a computer system for 
; predicting the pharmacokinetic behavior of a chemical substance in an insect, as well 
25 as a corresponding method and computer program product. 

The invention makes it possible, in a particularly efficient way, to 
calculate a prediction of the pharmacokinetic behavior of a chemical substance in an 
insect. In particular, the invention makes it possible to estimate the absorption, 
distribution and the excretion of chemical substances in insects on the basis of 
30 physicochemical parameters. 
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Summary of the invention 

To this end, a physiologically based phamiacokinetic simulation model 
of an insect is used to predict concentration/time profiles of a chemical substance in 
the compartments of the insect. The simulation model contains at least one parameter 
5 which is dependent on the substance to be studied. The parameter or parameters of 
the simulation model are, for a particular substance, predicted on the basis of one or 
more of the physicochemical properties of the substance. 

These physicochemical parameters are, for example, the lipophilicity of 
the substance, described by the distribution coefficient between water and 

10 phospholipid membranes or the octanol/water distribution coefficient, the molecular 
weight or the solubility. The relevant physicochemical parameters of the substance 
may either be determined by straightforward experiments or determined directly 
from this descriptor of the chemical structure of the substance by means of methods 
which are known per se, such as QSAR = Quantitative Structure Activity Relations, 

15 or neural networks. 

In the latter case, it is even possible to evaluate virtual substances by the 
method according to the invention, that is to say substances which have not yet been 
synthesized, with respect to their absorption properties and distribution properties in 
insects. Owing to the established relationship between physicochemical properties 

20 and pharmacokinetic properties, it is furthermore possible to derive general criteria 
for the optimization of insecticidal agents. 

According to a preferred embodiment of the invention, the rate 
coefficient of the iriter-compartmental mass transport, which is proportional to the 
product of the permeability for the substance and the effective svu-face area of the 

25 compartments, is used as a substance-dependent parameter of the simulation model. 
For each compartment of the simulation model, there is hence preferably a 
substance-dependent parameter which contains the permeability of the relevant 
compartment for the substance and the effective surface area of the relevant 
compartment. The permeability coefficient is a measure for the describing substance 

30 flux accross the cellular membrance. 

It is particularly advantageous that no further experimental studies are 
needed to determine the parameters for the simulation model based on a 
physicochemical property of the substance, but rather that the parameter or 
parameters can be determined on the basis of the physicochemical property of the 
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substance. This determination is carried out on the basis of a database, which has 
been previously determined experimentally for various test substances. The database 
contains the substance-dependent parameters of the simulation model which have 
been experimentally determined for the test substances, as well as the 
5 physicochemical properties of the test substances. This database is used to predict the 
substance-dependent parameter or paraineters for a substance to be studied. 

According to a preferred embodiment of the invention, a calculation 
function is obtained from the database by a linear regression. For example, the 
calculation function is a function of the lipophilicity and the molecular weight. To 
10 predict a parameter for the simulation model for a substance to be studied, it is hence 
merely necessary fpr the lipophilicity and the molecular weight of the substance to be 
studied to be evaluated with the calculation function, in order to obtain the 
parameter. With the predicted parameter, it is then possible to carry out a concrete 
simulation of concentration/time profiles during absorption and excretion of the 
15 substance in the insect. 

Alternatively, other prediction methods which are known per se may also 
be used instead of the calculation function obtained by linear regression. 

It is particularly advantageous that, after a database has been determined 
once for a particular number of test substances, no additional experiments are then 
20 needed for simulation of the pharmacokinetic behavior of further substances. This 
makes it possible to assess "candidates" for potential insecticides in respect of their 
pharmacokinetic behavior with a high throughput. Significant acceleration for 
research, development and optimization of new insecticides is therefore possible. 

Brief description of the drawings 
25 Preferred embodiments of the invention will be explained in more detail . 

below with reference to the drawings, in which: 

Figure 1 shows a block diagram of an embodiment of a computer system 
according to the invention. 

Figure 2 shows a flow chart of an embodiment of the method according 
30 to tihe invention, ' 

Figure 3 shows a physiologically based model of a caterpillar as a basis 
for a simulation model according to the invention with a substance-dependent 
parameter. 
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Figure 4 shows a diagram representing the database obtained for the 
simulation model of Figure 3 for carrying out a. linear regression, in order to 
determine a calculation function. 

Figure 1 shows a computer system 100. It may be a conventional 
personal computer (PC), a workstation or a client/server system. 

The computer system 100 has a simulation model 102. This is a 
physiologically based pharmacokinetic model of an insect. To this end, the 
simulation model l02 represents the compartments of the insect, and hence makes it 
possible to predict concentration/time profiles of a substance in the compartments of 
the insect. 

The simulation model 102 involves physiological parameters, 104 which 
depend only on the type of insect to be described, as well as one or more substance- 
dependent parameters 106. To carry out a simulation calculation of 
concentration/time profiles, it is hence necessary to input the parameter value of the 
substance-dependent parameter 106 for the substance to be evaluated. A concrete 
example of an embodiment of the simulation model 102 will be explained in more 
detail below with reference to Figures 3 and 4. 

The computer system 100 furthermore has a databank 108. The databank 
108 is used for storing a database, which have, been obtained on t6e basis of 
experimental studies of the pharmacokinetic behavior of test substances in the insect. 
For each of the test substances previously studied experimentally, the databank 108 
contains the parameter value or values which have been experimentally determined 
for the relevant test substance, as well as at least one physicochemical property of the 
relevant test substance. 

This database stored in the databank 108 forms the basis of the prediction 
of the parameter value for a new substance to be studied in the prediction module 
110. A calculation mle, which makes it possible to obtain the substance-dependent 
parameter value needed for carrying out the siniulation with the simulation model 
102 firom a physicochemical property of a substance to be studied, for example its 
lipophilicity or its molecular weight, is for example obtained by a linear regression 
method fi-om the database stored in the databank 108. 

The computer system 100 furthermore has an input/output module 112 
for input of the physicochemical property of a substance to be studied. The simulated 
concentration/time profiles are furthermore output via the input/output module 112. 
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The input/output module 112 may be coupled to a databank which 
contains descriptors of substances which actually exist or virtual substances, for 
example potential insecticides, and their physicochemical properties. In the case of 
virtual substances which have not yet been synthesized, the- physicochemical 
5 property required for input into the computer system 100 may be determined directly 
from the descriptor of the chemical structure of the test substance by means of 
methods which are known per se, such as QSAR, or neural networks. In the case of 
virtual substances, the descriptor instead of the physicochemical property may also 
be input via the input/output module 112. On the basis of the descriptor, the 
10 physicochemical property is then determined as an input quantity of the prediction 
module 1 10 in the computer system 100 itself. 

The simulation results are, for example, input into the database via the 
input/output module 112 so that evaluation of the simulation results can be carried 
out at a later tinie. 

15 Figure 2 shows a corresponding flow chart. In step 200, a 

physicochemical property of a substance to be studied is input. The physicochemical 
property is, for example, the lipophilicity or the molecular weight of the substance. 

In step 202, the parameter values of the substance-dependent^ parameters 
of the simulation model for the substance are calculated on the basis of a database 
20 based on the physicochemical property of the substance to be assessed. This database 
contains parameter values previously determined experimentally for various test 
substances. This calculation is carried out in step 202. 

In step 204, the parameters calculated in step 202 are input into the 
simulation model. The calculation of concentration/time profiles of the substance to 
25 be studied in the compartments of the insect is then carried out there. In step 206, 
these concentration/time profiles are output and can then be evaluated. 

In order to determine a database which is as meaningfiil as possible for 
prediction of the parameter values, it is advantageous if the test substances used to 
experimentally determine the database are as chemically diverse as possible. A new 
30 substance to be studied, for which the parameter values are to be predicted, should 
preferably lie within the test space represented by the test substances. 

Figure 3 shows, by way of example, a physiological model for a 
caterpillar. This physiological module is based on the model given by Greenwood et 
al. (see above). 
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The model involves the following eight compartments of the caterpillar: 
haemolymph (hi) 300, fat body (fb) 302, muscle (ms) 304, surface of the cuticle (cs) 
306, cuticle (c) 308, gut wall (gw) 310, gut content (gc) 312 and nerve bundle (nc). A 
substance to be studied is applied either topically via the cuticle 308, orally via the 
5 gut or by direct injection into tiie haemolymph 300. 

Like the blood in mammals, the haemolymph 300 is used as the main 
tr^sport phase between the various organs. 

Metabolism takes pace at the rate kc in the cuticle 308, at the rate khl in 
the haemolymph 300 and at the rate kgw in the gut wall 310. The excretion takes 
10 place at the rate kgc. 

The haemolymph 300 is modeled as a freely circulating liquid, which is 
in contact with the caterpillar organs described by the compartments. Transport 
processes between the compartments take place by passive diffusion, with 
permeation through the membranes as the rate-limiting step. The rate coefficients 
15 of the mass transport between the compartments are determined by the peraieability- 
surface area product PxAx and the volume Vx of the organs: 

■PA (1\ 
^x^-TT-^ (x€{c,mu,fb,nc,gw}) 

^X f 

f • 

where . 

Px = Permeability of the membrance (walls) of organ X for the substance, and 
20 Ax = surface area of the membrane (walls) of organ X. 

c = cuticle, 
niu = muscle, 
fb = fat body, 
nc = nerve bundle 
25 .gw = gut wall 

The equilibrium state is in this case reached after a time tx » Vx/(PxAx). 
The ratio of the concentrations in the peripheral compartments in the equilibrium 
state is in this case determined by the distribution coefficients related to the 
haemolymph Kx is defined as Kx^i: 

Kx=^ (xe{c,mu,fb,nc,gw}) * (2) 

30 " 
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Kx is hence the distribution coefficient of a substance between the 
haemolymph and the organ x in the equilibrium state. 

In addition, the distribution coefficients between the surface of the cuticle 
306 and the cuticle (Kc/cs) as well as the distribution coefficient between the gut wall 
310 and the gut content (Kgc/gw) must be known. 

The model with eight compartments which is known fi^om Greenwood et 
ai; has been modified in order to take account of 

- different types of application (oral, topical or by injection into the 
haemolymph), 

- an "open-loop" mass transport of food through the alimentary canal and 

- the time dependency of the organ volumes due to the larval growth. 

On the basis of this biophysical model of Figure 3, it is possible to 
describe a mass equilibrium rielation fpr each organ x of the caterpillar by a 
differential equation. Such mass equilibrium relations are given for the organs of the 
caterpillar in the appendix (Formulae Al to A8). 

In order to calculate a pharmacokinetic profile by means of the mass 
equilibrium relations of the organs of the caterpillar, a total of 23 parameter values 
must be known. They are 

- the organ volumes Vx (eight parameter values as well as additional 
parameter values for describing their time variations due to the larval growth), . 

- the distribution coefficients Kx of the substance between the 
haemolymph 300 and the peripheral compartments (c, mil, fb, nc and gw), 

- the distribution coefficients between the surface of the cuticle and the 
cuticle, as well as between the gut wall and the gut content, 

- the metabolic rate constants in the cuticle, the haemolymph and the gut 
wall, in the event that the substance to be studied is metabolized, and 

- rate coefficients Xx = [PxAx]A^x for the mass transport between the 
compartments (c, mu, fb, nc and gw). 

The distribution coefficients as well as the permeability-surface area 
products, and therefore the rate coefficients, depend on physiological parameters of 
the insect, as well as on the physicochemical properties of the substance to be 
studied. The metabolic rate constants are fiirthermore substance-specific. 
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The organ volumes can be experimentally determined by methods which 
are known per se. We can calculate the distribution coefficients Kx of the substance 
between the haemoljonph and each organ x as follows: 

• ^x/HjO 

~ with ^JitHfi - fmiur^ + ^/at ffai^ ^ pmin fpmeinjc (3) 

5 iri which 

Kfat = distribution coefficient of the substance in equilibrium between 
water and fat (lipophilicity) 

Kprotein ~ distribution coefficient of the substance in equilibriiun between 
water and protein 

10 fwater.x = volume fraction ofwater in compartment X, 

ffat,x = volume fraction of fat in compartment X, 

fprotein, x = volunie jfraction of protein in compartment X. 

The membrane affinity (MA) or, alternatively, the octanol/water 

distribution coefficient (Ko/w) of the substance can be used as an estimate of the 

15 distribution coefficients Kfat. The distribution coefficient Kprotein may, for example, be 

determined from the binding constants of human serum albumin (Kd"^"^ in [mmol]) 

and the molecular weight of the protein (65 kDa): ' , '* 

« 

Formula 3 can be simplified by combining the fat and protein fractions in 
20 a single organic fraction: 

^x/H O 

^x = — with ^x/H,0 = ^water.x + 0 " ^water,x) Kfat (5) 

^hl/HjO 

The rate constants of the mass transport between the compartments are 
experimentally determined for the test substances by matching them to experimental 
pharmacokinetic data and correlating them with the physicochemical properties of 
25 the test substances. 

As an approximation, it may be assumed that the permeability coefficient 
P for the substance is identical for all the organs x, in that is to say Px = P for all the 
organs x, with the exception of the cuticle. Furthermore, P is proportional to the 
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lipophilicity (Kfat = MA or Kfat = K<,/w, that is to say the octanol/water distribution 
coefficient) and the membrane diffusion coefficient Dmem of the substance: 

P<^Kf„D^ (6) 

The dependency of the membrane diffusion coefficient on the molecular 
5 . weight (MW) . of the substance is described by means-of an exponential relation: 

D^em ccMW-^-- . ; . (7) 

Formulae 6 and 7 can be combined with one another, so that 

P = a Kf3^ MW^'""- (8) 

results, in which a is a constant. By means of Equation 1, the rate constants of the 
10 ihass transport between the compartments can be expressed as follows: 

the following is obtained fi-om this by taking a logarithms 
Log{^) = Log(a^)-s_Log(MW) (lo) 

Figure 4 shows a corresponding logarithm/logarithm representation of 
15 Xoc/Kfat against the molecular weight MW of a plurality of test substances for yarious 
compartments. Figure 4 likewise shows the best-fit lines obtained by linear 
regressions through the measurement points. 

By means of linear regression, it is hence possible to determine the slope 
Smem as wcU as the intercept aAx/Vx from the experimentally determined database. A 
20 calculation rule has therefore be. obtained which makes it possible to calculate the 
quantity of a new substance, even a virtual one, merely on the basis of the 
lipophilicity (Kfat) and molecular weight MW. P* Ax is hence obtained on the basis of 
Formulae 8 and 9. On the basis of this substance-dependent parameter, as well as the 
equilibrium distribution coefficients Kx or Kc/cs and Kgc/gw from Eqs (3) and (5) 
25 respectively, it is hence possible to solve the system of equations in the appendix, for 
example numerically, and the time/concentration profiles of the substance to be 
studied in the individual compartments are obtained. 
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Appendix 

System of equations of the pharmacokinetic simulation model for the 

caterpillar 

(cs) ^cs^ = -fPAUC„-^) (Al) 

• ■ . . c/cs 

= (0 .^cr^ = [PA]e(:Ccs--^)+[PA]c(Cw-.^).-k,C^ - (A2)- 

c/cs . c 



(hi) 



dC C 

"d^^ = -ZfPAJx (Ch, --2-) - k^iCh, (X e {c.inu,fb.ncgw}) (A3) 



(n.u) .^'™'^ = t^^J-(^--r"> <A4) 



(fb) ^fl'^ = fPAJft.,(Chl-^) . ■ . (A5) 

(nc) '^-.^ = fP^J-(CH,-|^) . . (A6) 

dC C C 

(gc)- V^=[PA]gw(C«w--^) (A8) 

gc/gw . 



Cx: concentration of the substance in compartment x, 
: Vxi volume of compartment x; the volume may be time- variant according 
to an experimentally determined function (Vx(t)), 

[PA]x: permeability-surface area product of compartment X, 
10 Kx/y: distribution coefficient between compartments x.and y, 

Kx = Kx/hi: distribution coefficient between organ x and haemolymph, 
Kx: rate constants for the metabolism (x =c, hi, gw) 
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